div.frag=function(x,z)
{
pi=x[,2:ncol(x)]/apply(x[,2:ncol(x)],FUN=sum,MARGIN=1,na.rm=T)
H=-apply(pi*log(pi,base=exp(1)),FUN=sum,MARGIN=1,na.rm=T)

tabela.diversidade=data.frame(x[,1],H)

x11()
par(mfrow=c(2,3))

for(i in 1:ncol(z))
{
plot(tabela.diversidade$H~z[,i],xlab=names(z)[i],ylab="ndice de Shannon (H)",bty="l",cex=1,pch=16)
abline(lm(tabela.diversidade$H~z[,i]),col="red")
}
par(mfrow=c(1,1))
return(tabela.diversidade)
}



